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Introduction 

Solid  propellant  rocket  boosters  (SPRB)  are  used  extensively  in  aerospace  activities  when  a  strong 
thrust  is  needed  as  in  the  case  of  lifting  off  the  satellites  for  space  activities  [1],  Flow  structures 
occurring  within  the  SPRB’s  have  drawn  the  attention  of  the  researchers  either  experimentally  [2,3],  or 
as  it  happens  increasingly  today,  computationally  [4,5],  The  flow  field  inside  the  combustion  chamber 
experiences  different  processes  since  the  characteristic  velocity  throughout  the  chamber  is  on  the  order 
of  Ma  =  O(10"2)  [6],  whereas  the  exit  velocity  is  supersonic.  Time  dependent  behavior  of  various 
processes,  turbulence,  combustion,  as  well  as  with  the  compressibility  effects  add  intricacy  into  the 
simulation  processes.  Today,  major  study  field  for  the  SPRBs  is  the  flow  stability  problem  inside  the 
combustion  chamber  [6],  and  its  interaction  with  combustion  [7]  and  turbulence  [8], 

This  study  is  the  first  part  of  an  ongoing  research  held  in  the  Faculty  of  Mechanical  Engineering 
at  Istanbul  Technical  University  (ITU),  to  simulate  the  interaction  of  the  turbulence,  combustion  and 
the  acoustic  field  inside  SPRBs.  In  this  preliminary  work,  a  computational  study  has  been  performed  to 
identify  the  effect  of  various  turbulence  models  on  the  calculations  of  cold  flow  field  inside  a  model 
SPRB  combustion  chamber.  A  commercial  flow  solver,  Fluent  6.1.22,  has  been  employed  for  the 
computations. 

The  flow  configuration  chosen  is  an  idealization  of  that  found  in  a  solid  rocket  motor,  and  was 
selected  in  accordance  with  the  VECLA  facility  of  ONERA  that  is  an  experimental  set  up  for 
investigating  the  characteristics  of  injection  driven  flows  [4,8,9].  Schematic  of  the  facility  is  given  in 
Figure  1,  where  the  length  of  the  channel  is  L  =  0.581  m  and  its  height  is  h  =  0.0103  m.  The  channel  is 
bounded  at  y/h=  0  by  a  permeable  wall  allowing  the  mass  injection  and  at  y/h  =  l  by  an 
impermeable  wall.  The  upstream  at  x/L  =0,  head-end,  is  closed  and  the  downstream,  at  xl L  =  1,  exit 
section,  is  open  to  the  atmosphere.  Sidewall  mass  injection  is  used  to  mimic  the  normal  velocity  of 
gaseous  products  generated  by  combustion  of  gasified  propellant. 

For  this  particular  work,  based  on  the  Reynolds  averaged  Navier  Stokes  equations, 
computations  have  been  performed  by  modeling  the  turbulence  with  one  equation  model  (Spallart 
Allmaras),  two  equation  models  (Standard  k-e,  RNG  k-e.  Realizable  k-e,  SST  k-co )  and  the  Reynolds 
Stress  model.  Set  of  governing  equations  has  been  solved  by  employing  the  Finite  Volume  technique 
[10].  Boundary  conditions  were  specified  as,  air  injection  mass  flux,  m  =2.619  kg/m2s,  and  the 

appropriate  turbulence  quantities,  k  =0.0001  m  /s  ,  s  =0.001  m  /s  ,  co  =0.4  1/s,  at  the  permeable 
wall;  pressure  outlet,  137400  Pa,  in  the  exit  section  and  no  slip  condition  for  the  impermeable  wall  at 
y/h  =  \  and  the  solid  wall  in  the  head-end,  at  y/h  =  0.  Reynolds  number  based  on  the  injection 

velocity  and  the  channel  height,  Re; ,  is  approximately  8000.  Density  was  calculated  by  the  ideal  gas 

law,  and  hence  air  injection  temperature,  300  K,  was  specified.  Normalized  longitudinal  velocity  and 
the  turbulent  intensity  profiles  calculated  with  various  turbulence  closure  models  are  compared  with 
those  from  the  experiment  reported  in  [9]  at  several  cross  sections. 
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Figure  1.  Schematic  of  computational  domain 

Results  and  Discussion 

The  injection  Reynolds  number  was  chosen  as  to  produce  a  statistically  stationary  flow  field  and  the 
first  set  of  computations  dealt  with  the  comparison  of  the  longitudinal  velocity  profiles  obtained  by 
steady  and  time  marching  feature  of  the  FLUENT.  It  was  shown  that  the  profiles  match  excellently, 
allowing  us  to  exploit  the  steady  solver.  In  the  next  step,  several  turbulence  closure  models  have  been 
employed  for  computations  and  the  calculated  longitudinal  velocity  and  turbulence  intensity  profiles 
were  compared  with  those  obtained  from  an  experimental  study  reported  by  [9].  It  is  seen  from  Figure 
2  that,  Spalart-Allmaras  turbulence  closure  model  failed  to  compute  satisfactory  results.  Computations 
based  on  the  Reynolds  stress  model  predicted  the  normalized  longitudinal  velocity  profiles  in  a  good 
agreement  with  the  experiments.  However,  some  discrepancy  is  detected  in  the  neighborhood  of  the 
impermeable  wall,  hence;  the  validity  of  the  closure  model  is  questionable  in  this  region.  It  was 
possible  to  acquire  the  longitudinal  velocity  profiles  in  the  entire  channel  with  a  good  accuracy  with  the 
SST  k-co  model. 

As  given  in  Figure  3,  the  turbulence  intensity  levels  obtained  with  the  standard  k-e  model  failed 
to  capture  the  experimental  profiles.  However,  the  turbulence  intensity  profiles  obtained  with  the  SST 
k-co  model  are  satisfactory,  except  for  the  second  peak  observed  in  the  close  neighborhood  of  the  upper 
impermeable  wall.  Indeed,  the  accuracy  of  the  experiment  can  be  questionable  in  this  region  [9] 
because  the  existence  of  an  impermeable  wall  may  be  expected  to  increase  the  turbulence  levels  in  this 
region,  leading  to  a  possible  peak. 
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Figure  2.  Normalized  longitudinal  velocity  profiles,  at  x  =  0.500  m,  obtained  with  a)  Laminar  and 
Spalart-Allmaras,  b)  Standard  k-  s,  RNG  k-  e  and  Realizable  k-  s,  c)  Reynolds  stress  and  SST  k-co 

closure  models  and  from  the  experiment  f91. 
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Figure  3.  Turbulence  intensity  profiles  obtained  with  a)  Standard  k-  e,  b)  SST  k-co  and  from  c) 

Experiment  [9]. 
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Figure  1 .  Schematic  of  a  SPRB 


Figure  2.  Unsteady  flow  evolution  inside  a 
SPRB 
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Reynolds  Stress  Model 
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Computational  Domain 
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h 
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Permeable  Wall  -  Flow  Injection 


Schematics  of  the  Vecla  Facilty 


L  =  0.581  m  h  =  0.0103  m 


Common  Boundary  Conditions 
Injection  Surface:  Exit  Surface: 

m  =  2.619  kg/m2s  P  =  137400Pa 
T  =  300AT 

•Chaouat,  B.,  2000. 

•Kourta,  A.,  2004. 


Turbulence  Model  Specific 
Inlet  Conditions 

k-s  models  k  =  0.0001  m2 / s2 

s  =  0.00lm2/s 3 

k-co  model  k  =  0.000lm2 /s2 

co  =  0.4s”1 

RSM  /7  =  0.092  m2 /s2 
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Computational  Domain 
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TTT  T 


Finite  Volume  Method  (FLUENT) 
Simple  for  Pressure  velocity  coupling 
Second  order  upwind  for  remaining  equations 


100  x  290  mesh  points  in  x  and  y  for  steady  calculations 
100  x  242  mesh  points  in  x  and  y  for  unsteady  calculations 


International  Symposium  on  Energy  Conversion  Fundamentals,  21-25  June  2004,  Istanbul  i6 


Longitudinal  Velocity  [m/s]  at  x=0. 15  m 


Time  Marching  Approach 
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Is  1.0000 
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0.030  0.060  0.090  0.120  0.150  0.180  0.210 

Time  [sec] 


Time  history  for  longitudinal  velocity  at  sections  a)  0.150  m, 
b)  0.55  m  at  mid  height,  and  c)  outlet  to  inlet  mass  flux  ratio,  obtained  by  time 

marching  approach. 


SST  k-to  model  has  been  used  for  time  marching  computations 
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Channel  Height  [m] 


Time  Marching-Steady  Approaches FacultyofMechanical 
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Normalized  Longitudinal  Velocity 


a) 


b) 


c) 


d) 


Comparison  of  the  normalized  longitudinal  velocity  distributions  obtained  by  time  marching  and 
steady  schemes  at  sections  a)  0.031  m,  b)  0.120  m.,  c)  0.220  m.,  d)  0.350  m.,  e)0.400  m.,  f)  0.450  m., 

g)  0.500  m.  and  h)  0.570  m.  for  the  coarse  mesh. 


SST  k-to  model  has  been  used  for  comparison 
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Channel  Height  [m] 


Laminar  and  Spalart  Allmaras 
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Normalized  Longitudinal  Velocity 


a) 


b)  c)  d) 


Normalized  longitudinal  velocity  profiles  obtained  with  the  Laminar  approach  ,  the  Spalart 
Allmaras  turbulence  closure  model  and  from  the  experiment  [3]  at  sections  a)  0.031  m,  b)  0.120 
m.,  c)  0.220  m.,  d)  0.350  m.,  e)0.400  m.,  f)  0.450  m.,  g)  0.500  m.  and  h)  0.570  m.  for  the  fine  mesh. 
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Computations  with  k-z  based 
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Normalized  Longitudinal  Velocity 


a) 


b)  c)  d) 


Normalized  longitudinal  velocity  profiles  obtained  by  the  Standard  k-  z,  RNG  k-z,  Realizable  k-z 
turbulence  closure  models  and  from  experiment  at  sections  a)  0.031  m,  b)  0.120  m., 
c)  0.220  m.,  d)  0.350  m.,  e)0.400  m.,  f)  0.450  m.,  g)  0.500  m.  and  h)  0.570  m.  for  the  fine  mesh. 
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Reynolds  Stress  Model 
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Normalized  longitudinal  velocity  profiles  obtained  with  the  Reynolds  Stress  turbulence  closure 
model  and  from  experiment  at  sections  a)  0.031  m,  b)  0.120  m.,  c)  0.220  m.,  d)  0.350  m.,  e)0.400  m., 

f)  0.450  m.,  g)  0.500  m.  and  h)  0.570  m.  for  the  fine  mesh. 
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Channel  Height  [m] 


SST  k-omega 
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Normalized  Longitudinal  Velocity 


a) 


b) 


c)  d) 


Normalized  longitudinal  velocity  profiles  obtained  by  the  SST  k- w  turbulence  closure  model  and 
from  experiment  at  sections  a)  0.031  m,  b)  0.120  m.,  c)  0.220  m.,  d)  0.350  m.,  e)0.400  m.,  f)  0.450 

m.,  g)  0.500  m.  and  h)  0.570  m.  for  the  fine  mesh. 
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Turbulence  Intensity 
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Turbulence  intensity  profiles  obtained  with  the  a)  Standard  k-e,  b)  RNG  k-e,  c)  Realizable  k-e, 
d)  Reynolds  stress,  e)  SST  k-cj  turbulence  closure  models  and  f)  from  experiment,  at  eight  sections 
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Results 
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•  k-e  based  computations  have  not  produced  velocity  and  turbulence 
intensity  profiles  compatible  with  the  experimental  profiles,  neither  in 
the  pre  (laminar)  nor  post-transition  (turbulent)  sections. 

•  Computations  based  on  the  Reynolds  Stress  Model  have  predicted 
the  normalized  longitudinal  velocity  profiles  well  except  for  the 
laminar  region,  and  turbulence  intensity  profiles  are  similar 
qualitatively  as  those  from  the  experiments. 

•  SST  k-co  turbulence  closure  model  has  yielded  the  longitudinal 
velocity  profiles  in  a  good  agreement  with  the  experiment  [3]. 
Furthermore,  accuracy  of  the  turbulence  intensity  profiles  was 
satisfactory  compared  with  the  results  from  the  rest  of  the  closure 
models. 
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